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ABSTRACT 

The total infrared (TIR) luminosity from galaxies can be used to examine both star formation and dust 
physics. We provide here new relations to estimate the TIR luminosity from various Spitzer bands, in particular 
from the 8 pm and 24 pm bands. To do so, we use 45" subregions within a subsample of nearby face-on spiral 
galaxies from the Spitzer Infrared Nearby Galaxies Survey (SINGS) that have known oxy gen abundances as 
well a s integrated galaxy data from the SINGS, the Local Volume Legacy Survey (LVL) and lEngelbracht et alj 
( 120081) samples. Taking into account the oxygen abundances of the subregions, the star formation rate intensity, 
and the relative emission of the polycyclic aromatic hydrocarbons at 8 pm, the warm dust at 24 pm and the 
cold dust at 70 /im and 160 /im we derive new relations to estimate the TIR luminosity from just one or two of 
the Spitzer bands. We also show that the metallicity and the star formation intensity must be taken into account 
when estimating the TIR luminosity from two wave bands, especially when data longward of 24 /im are not 
available. 

Subject headings: infrared: galaxies 



1. INTRODUCTION 

With the end of the Spitzer cold phase and the widespread 
availability of 8 pm and 24 /im bands observations in the 
archives, the availability of relations to determine the total in- 
frared (TIR) emission from these wave bands by themselves 
is crucial to efficiently exploit the archives. In addition, the 
Herschel Space Observatory observes dust emission at rest- 
frame 24 /im and longer wavelengths for galaxies redshifted 
to z — 1.5 in the PACS 60 /im and longward bands with a 
spatial resolution as good as the Spitzer 24 pm band, for in- 
stance. Thus, estimating the total infrared flux from these 
bands becomes crucial for measuring total infrared fluxes us- 
ing Herschel data. Finally, the advent of new instrumentation 
in the coming years, such as the James Webb Space Telescope 
(JWST) or the Atacama Large Millimeter Array (ALMA) will 
also open a new window on the infrared emission of nearby 
and distant galaxies. For instance, the 18 /im JWST/MIRI 
band and the 350 /im ALMA band will probe the rest-frame 
8 /im and the 160 /im emission of z — 1.2 galaxies. 

Determining the total infrared emi ssion using the 8-160 /im 
Spitzer bands wit h equation 4 from Dale & Heloul d2002l) or 
equation 22 from lDraine & Lil (120071) yields a better estimate 
of the total infrared flux than using a single wave band as a 
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proxy for the total infrared flux, but they necessitate 3 and 4 
infrared bands respectively. Indeed, these relations necessi- 
tate using the much lower resolution 70 /im and 160 pm data, 
which have resolutions of 18" and 40", respectively. The 
poorer resolution of the 70 pm and 160 pm bands strongly 
constrains the scales on which total infrared fluxes can be 
measured, even for local galaxies. Attempts to derive a re- 
lation to estimate the total infrared em ission from th e 8 /im 
and 24 pm bands have been made by ICalzetti et al.l (12005 ) 
using NGC 5194 (M51) data, iPerez -Gonzalez et alj (12006ft 
using NGC 3031 (M81) data, and iThilker et al.l (120071) using 
M33 data. For each galaxy, th e scatter around the relation is 
about 40%. iRieke et ail (120091) recently showed that using the 
24 pm band only provided good re s ults. U sing the 24 pm and 
the 70 pm bands IPapovich & Belli (120021) have also provided 
an estimate of the TIR emission with a similar uncertainty. 
However, these relations may be applicable only to galaxies 
with similar metallicities and star formation rate intensities 
(Calzetti et al. 2010, submitted; Li et al. 2010, in prepara- 
tion). Metallicity variations have been associated with varia- 
tions in mid- and far-infr ared colors, as has been observ ed in 
metal-po or galaxies (e . g. lEngelbracht et al.ll2005l 120081) . In- 
deed, the ICalzetti et al.l (120051) relation, which was derived us- 
ing data from a very metal-rich galaxy, underestimates the to- 
tal infrar ed emission by a factor of a few in metal-poor dwarf 
galaxies dCannon et 
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In this article we derive relations to estimate the total in- 
frared emission using just one or two of the Spitzer bands, 
with a strong focus on deriving the total infrared flux from 
just the higher resolution, shorter wavelength 8 pm and 24 pm 
bands. To do so, we use regions within a subset of nearby 
face-on spiral galaxies from the Spit zer Nearby Galaxies Sur- 
vey (SINGS, iKennicutt et all 120031) that are resolved in the 
Spitzer 160 pm band as well as integrated galaxy luminosities 
for galaxies in the SINGS, the Local V olume Legacy survey 
(LVL) and the lEngelbracht et al.l (120081) . hereafter E08, sam- 
ples. In section|2] we describe the samples of galaxies and the 
data processing. In section @] we present the results, and we 
discuss them in section [5] Finally we summarize our results 
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and conclude in section [6] 

2. SAMPLE AND DATA 

We use several samples to derive the TIR emission both 
from galaxy subregions and g alaxies. We use t he subset 
of SINGS galaxies selected by iBendo et al] (120081) to study 
galaxy subregions. The requirements to perform this study 
are similar to theirs: the galaxies must have major axes of 
5' or greater so that galaxy substructure is resolved; the in- 
clination must be no more than 60° from face-on or less so 
that spatial variations can be seen with little overlap due to 
projection effects; and spatially resolved oxygen abundance 
data must be available. Galaxies which have compact infrared 
emission (NGC 1512, NGC 4826 [M64]), where muxbleed 
effects at 8.0 pm caused problems with interpreting the data 
(NGC 1097, NGC 1566 and NGC 4736 [M94]), or where very 
bright foreground stars caused problems with interpreting the 
IRAC data (NGC 3561) were exc luded. NGC 393 8 and 
NGC 4579 (M58), which were in the lBendo et all (120081) sam- 
ple, were not used here. The optical spectrum of NGC 3938 
is very noisy, which makes the oxygen abundance measure- 
ments unreliable, while the optical spectrum of NGC 4579 
is strongly affected by an AGN that causes problems when 
determining the oxygen abundances (Moustakas et al. 2010, 
in preparation). The final sample consists of 13 galaxies: 
NGC 0628 (M74), NGC 0925, NGC 2403, NGC 3031 (M81), 
NGC 3184, NGC 3351 (M95), NGC 4254 (M99), NGC 4321 
(M100), NGC 4725, NGC 5055 (M63), NGC 5194 (M51), 
NGC 6946 and NGC 7793. 

Observations and da ta processing information are provided 
by iBendoe t al. (2008). The extraction of flux densities for 
subregions in these galaxies was performed using the method 
described by IBendo et alj (12008b . First, the data were c on- 
volved with kernels developed by iGordon et al.l (12008b to 
match their PSF to that of the 160 /im data. The entire analysis 
is performed on the data smoothed to the 160 /im image res- 
olution. Next the stellar continuum emission was subtracted 
from the 8.0 /im and 24 /im bands using the relations 

L v (PAH 8 pm) = L v (8 pm)- 0.232 x L„(3.6 pm) (1) 

L v (24 pm) = L v (24 pm obs)-0.032 x L J/ (3.6 pm) (2) 



given bv lHelou et ail d2004l) . The images were then rebinned 
into 45" pixels, and pixels with low S/N and pixels strongly 
affected by artifacts or foreground stars were masked out. The 
processed images of NGC 3031 are shown in Figure Q] as an 
example. The flux densities in these 45 " pixels are then used 
in the analysis. See IBendo et al.l (120081) for additional details. 
The physical resolution ranges from 0.7 kpc to 3.6 kpc with 
a mean of 2.0 ±1.0 kpc. The oxygen abundances have been 
calculated for each 45" region using the abundance gradients 
from Moustakas et al. (2010, in preparation), assuming that 
the gradient in each galaxy is azimuthally symmetric. They 
range from 8.29 to 8.93 with (1 2 + logO/i/) = 8.60±0 . 13 av- 
eraging the estimates fro m the |Kobulnickv & Kewlevl (120041) 
(hereafter KK04) and the lPilvugin & Thuanl (120051) (hereafter 
PT05) estimators. Indeed these estimators are representative 
of the maximum and minimum oxygen abundance and as such 
bracket the actual ab undance. A similar m ethod was used for 
the same reasons by iCalzetti et alj ( 120071) . Unless specified 
otherwise we use this average throughout the paper to esti- 
mate the metallicity of galaxy subregions. 

We use integra ted galaxy data f r om the SINGS (|Dale et al.l 
120071) and LVL dLee et al.l 120081; iDale et all 120091) samples 
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FIG. 1 .— Processed PAH 8, 24, 70, and 160 fim images of NGC 3031. The 
intensity of the color is proportional to the flux. Each pixel has an angular 
resolution of 45". North is up and east is left. 

as well as the E08 sample as three additional data sets in 
this analysis. Although only half of the LVL galaxies cur- 
rently have oxygen abundances available (Marble et al. 2010, 
submitted), these galaxies along with the SINGS and E08 
galaxies allow us to probe a larger range of metallicities 
than the SINGS galaxies subregions. The SINGS sample 
contains galaxies with a metallicity from 8.02 to 8.99 with 
(12 + logO /j/) = 8.64 ± 0.21 a veraging the KK04 and PT05 
estimators (Cal zetti et al.ll2007h . Finally, the E08 sample "in- 
cludes well-known starbursting or star-forming galaxies from 
the literature". It probes a large range of oxygen abundances, 
from 7.31 to 8.85, using electron t emperature measurements 
and 03N2 dPettini & Pagell 12004b for the most metal-rich 
galaxies. The sampling is not even and is mainly concen- 
trated on_Jow^er : mid_^uridances: (12 + logO/H) = 8.28 ± 
0.33 jEngelbracht et al.11200 8). Galaxies which were not de- 
tected in one or more bands are excluded from this study. In 
addition, all galaxies are selected with S/N>5. The final sam- 
ple contains 57 of the 75 SINGS galaxies, 179 of the 258 LVL 
galaxies and 48 of the 66 E08 galaxies. 

3. METHODS TO DERIVE THE RELATIONS 

As a first step, we estimate the total infrared luminosity 
from the PAH 8 nm to 160 pm Spitzer bands. To do so we 
use equation 22 fromlDraine &11 d2007l) - a newer derivation 
than IDale & Heloul (120021) - based on models that have been 
calibrated using results from Spitzer. 



L(TIR) = 0.95L(PAH &pm)+1.15L(24 pm) 
+L(70 pm)+L(l60 pm), 



(3) 



where L = v h v . The difference with equation 4 of 



IDale & Heloul (12002b is minimal, 0.03 ±0.03 dex (Figure 
and section|5]l. 

To derive the relations to determine the total infrared lu- 
minosity from the combination of Spitzer bands we proceed 
in two ways. The first way is to perform a linear fit on the 
colors. For instance we calculate the coefficients a and b 
of a relation of the form: log (L(TIR)/L(24 pm)) = a + b x 



log (L(PAH 8 [mi) / L(24 /im)) . For easier use, we provide the 
relations under the form: logL(TIR) = logL(24 [im) + a + b x 
log (L(PAH 8 jttm)/L(24 fjm)). We use distance independent 
quantities - the ratio of two luminosities - in order to avoid 
the well known correlation bias induced by luminosity ver- 
sus luminosity relations. Indeed luminosity versus luminosity 
correlations are observed whether or not the data are actu- 
ally correlated due to the multiplication of the flux density 
by the square of the distance. The drawback of such a rela- 
tion is that the luminosities in the two Spitzer bands are not 
taken into account independently but are tied through the pa- 
rameter b. To alleviate this limitation, we have performed 
fits in another distance independent way using the luminos- 
ity per unit area S. The relations obtained are of the form: 
log£(ra?) = a + b x logT,(PAH 8 fim) + c x log £(24 /zra) for 
instance. One drawback of this formulation is that it requires 
the target galaxy to be resolved - which is often not the case 
for distant galaxies, especially in mid- and far- infrared bands 
- and to know its distance. 

In addition to deriving the relations, we also study the in- 
fluence of the metallicity and the star formation intensity (by 
way of the total infrared luminosity per unit area). To do so 
we perform the aforementioned fit dividing the sample into 
5 bins each containing the same number of data points in 
12 + logO/H or £(7Yr?). Evaluating the slope change as a 
function of the abundance allows us to study in detail its in- 
fluence on the determination of the TIR luminosity. Finally, 
we also provide relations of the form log£(7Yr?) = a + b x 
logY,(PAH 8 /im) + c x log £(24 fj,m) + d x (l2 + \ogO/H) to 
correct for the effect of the metallicity on the estimation of the 
total infrared luminosity per unit area. 

4. RESULTS 

The relations obtained through the methods described in the 
previous section are listed in Tables[T|[2]and[3] We present the 
results hereafter. 

4.1. TIR estimates from the combination of the 8 and 24 fim 

bands 

As PAH and hot dust are significant contri butors to the 
total infrared luminosity in metal-rich galaxies dDraine & Lil 
l2007h . the 8 fim and 24 fim observations permit an accurate 
determination of the TIR luminosities. Furthermore, as their 
spatial resolution is significantly better than the one of the 
70 fim and the 160 fim bands, resolving local galaxies is pos- 
sible. 

4. 1 . 1 . Metallicity effects 

It is well known that the presence of dust carriers influ- 
encing the 8 /im band emission is directly affected by the 
metallicity (lEngelbracht et al.l 12.005b iRosenberg et al.l 1 2006 ; 
WuetalJ 120061: : [Madden et all 12006b I Jackson et all 12006 : 
Draine & Lil2007l) . The large number of galaxy subregions in 
the sample studied here allows us to take into account accu- 
rately this parameter in the determination of the total infrared 
emission from the PAH 8 /im and 24 /im bands. In addition, 
the integrated SINGS and the E08 galaxies have oxygen abun- 
dances that are averaged over all the galaxy. 

The relations between the L(24 fim)/L(TIR) and 
L(PAH 8/im)/L(24 /im) colors for all the samples are 
presented in Figure [2] In order to clearly show the difference 
in the behavior of integrated galaxies and galaxy subregions, 
we plot them together. For comparison we also show the 
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rel ations for NG C 5194 , NGC 3031, and M33 determ ined 
bv ICalzetti et al.1 (12005b . iPerez-Gonzalez et al.l (120061) . and 
iThilker et all (120071) res pectively. The NGC 5194 relation 
by ICalzetti et al.l (120051) is offset from the other relations 
and also from the observations as it has a much higher 
L(24 fim)/L(TIR) ratio. The reason is that the NGC 5194 
flux densities were measured after removing a local back- 
ground, whereas for NGC 303 1 and other data, flux densities 
were measured after removing only a global background. 
The difference therefore could be an indication of the relative 
contribution of the diffuse large scale infrared emission 
compared to point-like sources. The parameters of the best 
fitting relations are given in Table [T] separately for each 
sample and the galaxy subregions. In addition, we provide 
relations between the total infrared luminosity per unit area 
as a combination of the luminosity per unit area in the PAH 
8 /jm and 24 fim bands in Table [2] and we also take into 
account the oxygen abundance in Table [3] 




-1.0 -0.5 0.0 0.5 

log(L(PAH 8um)/L(24)) 



FIG. 2. — Ratio of the 24 fim to the total infrared emission - as calcu- 
lated using equation [5]- versus the ratio of the PAH 8 fim to 24 fim emis- 
sion. The "x" symbol represents the data for the subregions in the sam- 
ple of face-on galaxies. The color of the "x" symbols represent the oxygen 
abundances as given by Moustakas et al. (2010, in preparation). The filled 
circles represent integrated galaxy data for the SINGS sample, the empty 
circles represent integrated galaxy data for the LVL sample, and the filled 
triangles represent integrated galaxy data for the E08 sample. The solid 
lines represen t various relat i ons: t he yellow one for the relation found in 
NGC 5194 by Calzetti et al. (2 005T). the mag enta one for the relation found 
in NGC 3031 by Perez-Gonzalez et al. 12006), the purple one for the relation 
found in M 33 by 'Thilker et al. 12007), the blue one for LVL galaxies, the 
red one for SINGS galaxies, the green one for E08 galaxies, and the black 
one for galaxy subregions. For the relations published in the literature, the 
L(PAH 8 fim)/L(24 fim) range over which they were derived are displayed 
with a solid line, and a dashed line out of these bounds (P. Perez-Gonzalez 
and D. Thilker, private communications). 

First of all in Figure [2] we notice that despite the scatter, 
different relations can be seen for subregions with different 
abundances. Indeed, we see that the slope is clearly shallower 
for the lower oxygen abundance subregions which describe a 
lower branch in the diagram. Interestingly, the low metallic- 
ity trend is followed by a part of the SINGS and most of the 
LVL samples. Even though the LVL sample spans a signifi- 
cant range of metallicities, it is chiefly constituted of dwarf 
galaxies. As a consequence, LVL galaxies are statistically 
more metal-poor than SINGS ones and therefore exhibit for 
a significant number of them little or no PAH emission. So in 
this particular case the data follow a trend similar to the one 
for the lowest oxygen abundance galaxy subregions (note that 
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TABLE 1 

TIR ESTIMATIONS FROM LUMINOSITIES 



1.012 ±0.009 
0.997 ±0.023 
1.141 ±0.010 
0.887 ±0.034 
4.580±0.108 
6.696 ±0.532 
8.519 ±0.504 
7.075 ±0.790 
4.961 ±0.109 
3.810±0.975 
1.838±0.326 
3.487 ±1.022 



0.510±0.021 
0.494 ±0.071 
0.120 ±0.024 
0.400 ± 0.074 
0.888 ±0.003 
0.836 ±0.015 
0.782±0.015 
0.829 ±0.022 
0.887 ±0.003 
0.923 ±0.028 
0.979 ±0.010 
0.924 ±0.029 



0.096 
0.142 
0.132 
0.185 
0.080 
0.113 
0.218 
0.182 
0.082 
0.194 
0.157 
0.212 



y *i X2 

L(TIR SU b) L(24fim) L(PAH 8fi?n) /L(24 fim) 

L(TIR sin gs) L(24fim) L(PAH %fim) /L(24 fim) 

L(TIR LVL ) L(24fim) L(PAH Sfim) /L(24 fim) 

L(TIR E0& ) L(24fim) L(PAH 8 fim) /L(24 fim) 

L(TIR SV b) ■■■ L(PAH 8 fim) 

L(TIRsings) ■■■ L (PAH 8 fim) 

L (TIRlvl) ■■■ L (PAH &fim) 

L (TIR E0S ) ■■■ L (PAH %fim) 

UTIRsub) ■■■ L(24iim) 

WIRsings) ■■■ L(24fim) 

L(TIR LVL ) ■■■ L(24fim) 

L(TIR EOg ) WAjjm) 

L(TIR S ub) L(70 fim) L (PAH 8 fim) /L (70 fim) 

L(TIR SING s) WO /im) L(PAH 8 fim) /L(10 fim) 

L (TIR LVL ) L (70 fim) L (PAH 8 fim) /L (70 fan) 

L(TIR E0S ) L(l() fim) L(PAH 8 fim) /L(70 fim) 

L(TIR SUB ) L (70 fim) L (24 fim) /L (70 fim) 

L(TIR SINGS ) L(10fim) L(24 fim) /L(10 fim) 

L(TIR LV l) WO fim) L (24 fim) /L (70 fim) 

L (TIR E0S ) L (70 fim) L (24 fim) /L (70 fim) 

WIRsub) ■■■ WO fim) 

WIRsings) ■■■ WO fim) 

WIRlvl) ■■■ L(10fim) 

L(TIR Em ) L (70 fim) 

L(TIR SU b) L(\60fim) L(PAH 8 fim) /L(160 fim) 0.454 ±0.005 0.283 ±0.010 0.050 

L(TIR SINGS ) L(\60fim) L(PAH 8 fim) /L(160 fim) 0.599±0.037 0.370±0.069 0.125 

L(TIRlvl) L(\60fim) L(PAH 8 fim) /L(160 fim) 0.458±0.019 0.027±0.023 0.118 

L(TIR EOi ) L(\60fim) L (PAH 8 fim) jh (160 fim) 0.831 ±0.038 0.310±0.083 0.182 

WIRsub) L (160 fim) L (24 fim) /L (160 fim) 0.616 ±0.005 0.332 ±0.006 0.031 

UTIRsings) L (160 fim) L(24 fim)/L( 160 fim) 0.714±0.014 0.436±0.019 0.047 

L(TIR LV l) L(l60fim) L (24 fim) /L (160 fim) 0.751 ±0.013 0.447 ±0.017 0.055 

L(TIR E0S ) L(l60fim) L(24 fim)/L( 160 fim) 0.748 ± 0.010 0.465 ± 0.024 0.070 

L(TIR SUB ) L(l60fim) L (70 fim) /L (160 fim) 0.408 ±0.002 0.390 ±0.006 0.030 

WIRsings) L(l60fim) L(10fim)/L(\e0fim) 0.437 ±0.007 0.623 ±0.030 0.052 

L(TIR LV l) L(160fim) L (70 fim) /L (160 fim) 0.395 ±0.003 0.552 ±0.015 0.040 

L(TIR E0S ) L(l60fim) L (70 fim) /L (160 fim) 0.393 ±0.025 0.854±0.054 0.054 

L(TIR SU b) ■■■ L(l60fim) -1.278±0.093 1.047±0.003 0.057 

L(TIR S ings) ■■■ L(l60fim) 1.683 ±0.851 0.965±0.024 0.159 

L(TIR L vl) ■■■ L(l60fim) 1.642±0.263 0.965±0.008 0.114 

L(TIR EOg ) ■■■ L(\60fim) 1.747± 1.004 0.971 ±0.028 0.215 

NOTE. — Coefficient for the fit logy = logxi +a + b X logJC2- <r is the standard deviation of the data 

points around the best fit line. The subscripts refer to the sample. 



0.662 ± 0.003 
0.583 ±0.017 
0.486 ±0.010 
0.478 ±0.018 
0.789 ±0.013 
0.497 ± 0.050 
0.583 ±0.034 
0.415 ±0.020 
1.210±0.142 
-1.057 ±0.580 
-1.210±0.222 
-1.098 ±0.409 



0.423 ±0.010 
0.300 ±0.032 
0.158=1=0.011 
0.198 ±0.023 
0.351 ±0.020 
0.078 ±0.075 
0.287 ±0.043 
0.183 ±0.039 
0.981 ±0.004 
1.042 ±0.016 
1.045 ±0.006 
1.040 ±0.011 



0.059 
0.063 
0.067 
0.055 
0.083 
0.101 
0.089 
0.073 
0.095 
0.101 
0.076 
0.078 



12+logO/H > 8.29 for the subregions). For most of the LVL 
galaxies, the L(PAH 8/im)/L(24 fim) ratio is primarily depen- 
dent on the PAH 8 /im luminosity which itself is dependent on 
the metallicity. That is, the L(24 nm)/L(TIR) ratio is hardly 
affected by a change of the L(PAH 8 /im) j L(24 fim) color. 
Quantitatively, we see in Table[T]that the slope is indeed much 
shallower for SINGS and especially LVL galaxies, which 
show little dependence on L{PAH 8 fim) /L(24 fim) as could be 
expected. The subregions sample has the highest mean slope 
due to being constituted of spiral galaxies only. The behavior 
of the E08 sample shows that the metallicity is not the only 
parameter driving the correlation between L(24 fim)/L(TIR) 
and L(PAH %fim)/L{2A fim). Most of the sample seems to 
follow the trend set by higher metallicity subregions describ- 
ing the upper branch in the diagram. As this sample has been 
specifically constituted with star-forming and starburst galax- 
ies, the chief parameter that drives the variation could be the 
intens ity of star formation. We will probe its effect in section 
ECU 

In Table |3] we estimate the total infrared luminosity per unit 
area as a function of the luminosity per unit area in the PAH 
8 fim and 24 fim bands as well as the oxygen abundance. We 
see that combining the PAH 8 fim and 24 fim bands, the con- 
tribution of the metallicity is smaller than when estimating 



from the PAH 8 fim band only. The reason is that the combi- 
nation of these two bands is an indirect measure of the metal- 
licity which permits to obtain a more accurate measurement 
than when no information on the metallicity is available. In- 
deed, taking into account the oxygen abundance, averaging 
over the subregions, SINGS and E08 samples, the scatter is 
reduced by 0.003 dex when combining the PAH 8 fim and 
24 fim bands but 0.025 dex when estimating the total infrared 
emission from the one in the PAH 8 fim band only. 

In order to study how the fit parameters evolve as a func- 
tion of the oxygen abundance, in Figure [3] we plot the best fit 
parameters for different bins of oxygen abundances following 
the method described earlier in section [3] 

Three different trends can be seen with increasing oxygen 
abundance: 1. at low oxygen abundance, the slope is shallow, 
2. at a higher oxygen abundance, the slope is steeper, this 
shows the greater dependence on the PAH emission, 3. finally 
at the highest oxygen abundance, the slope decreases slightly. 
Considering that L(24 fan) /L(TIR) is a f unction of the in- 
tensity of the illu minating radiation field dDale et al.l 12.001b 
iDraine & Lill2007l) . the obse rved trend reflects the results ob- 
tained by oth er authors (e. g. lCalzetti et a l. 2005; Bendo etlril 
120061 120081) that the relative strength of PAH emission de- 
creases as the strength of the illuminating radiation field in- 
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TABLE 2 

TIR ESTIMATIONS FROM LUMINOSITIES PER UNIT AREA 



T,(TIR SUB ) 

^(TIRsings) 
Z(TIR LVL ) 
£(77P £08 ) 
Z(TIR SUB ) 

£ (TIRsings) 
T,{TIR LVL ) 
X(TIR E0S ) 
Z(TIR SUB ) 

Z (TIRsings) 
T,(TIR LVL ) 
Z(TIReoss) 
V(T1R S ub) 

£ (TIRsings) 
Z(TIR LVL ) 
£(77P £ os) 
£(77/W 

£ (TIRsings) 
Z(TIRlvl) 

£(77R £0 8) 
X(TIR SUB ) 

T,(TIR s ,ngs) 

£(77Klvl) 

£(77R £0 8) 



£ (PAH 8/im) 
£ (PAH 8/im) 
£ (PA// 8/im) 
£ (PA// 8/im) 
£ (PA// 8fim) 
£ (PA// 8fun) 
£ (PA// 8(U/n) 
£ (PA// 8fim) 
£(24 Aim) 
£(24 ftm) 
£(24 fim) 
£(24 fim) 
Y.(PAH'ifim) 
£ (PA// 8 fin i) 
£ (PA// 8 fim) 
£ (PA// 8 fim) 
£ (24 fim) 
£ (24 ftm) 
£(24 (im) 
£ (24 fim) 
£ (70 fim) 
£ (70 fim) 
£ (70 fim) 
£ (70 fim) 



*2 

£ (24 fim) 
£ (24 fim) 
£ (24 ftm) 
£ (24 fim) 



£(70 /iin) 
£(70 ftm) 
£ (70 fim) 
£(70 /im) 
£ (70 fim) 
£ (70 fim) 
£ (70 fim) 
£ (70 fim) 



5.692 ±0.124 
2.085 ±0.380 
5.408 ±0.416 
6.795 ±0.780 
5.982 ±0.170 
2.766 ±0.476 
13.085 ±0.663 
7.771 ±1.515 
6.418 ±0.169 
1.714 ±0.540 
4.220 ±0.466 
8.083 ±0.898 
3.994 ±0.070 
0.627 ±0.207 
1.978 ±0.295 
0.085 ±0.433 
4.593 ±0.137 
-0.018±0.319 
1.078 ±0.400 
1.260 ±0.698 
3.847 ±0.176 
-0.1 18 ±0.299 
0.042 ±0.417 
-0.617 ±0.673 



0.433 ±0.014 
0.525 ±0.067 
0.174 ±0.020 
0.257 ±0.053 
0.845 ±0.005 
0.944 ±0.014 
0.629 ±0.021 
0.799 ±0.045 
0.842 ±0.005 
0.981 ±0.016 
0.904 ±0.014 
0.783 ±0.027 
0.425 ±0.006 
0.302 ±0.034 
0.178±0.011 
0.194 ±0.023 
0.405 ±0.015 
0.069 ±0.076 
0.3 10 ±0.045 
0.216 ±0.047 
0.902 ±0.005 
1.017 ±0.009 
1.010±0.013 
1.028 ±0.020 



0.425 ±0.014 
0.442 ± 0.070 
0.693 ±0.026 
0.566 ±0.050 



0.476 ±0.006 
0.697 ±0.037 
0.777 ±0.016 
0.817 ±0.028 
0.483 ±0.016 
0.946 ±0.078 
0.675 ± 0.049 
0.760 ±0.061 



0.062 
0.133 
0.105 
0.122 
0.086 
0.175 
0.233 
0.242 
0.087 
0.192 
0.127 
0.151 
0.033 
0.063 
0.063 
0.054 
0.063 
0.100 
0.089 
0.072 
0.082 
0.101 
0.100 
0.087 



£(77/W 


£ (PAH 8 fim) 


£(160 


fim) 


-1.794±0.180 


0.156±0.014 


0.908 ±0.018 


0.046 


X(TIR s ,ngs) 


£ (PAH 8 fim) 


£(160 


fim) 


1.474 ±0.365 


0.408 ± 0.068 


0.567 ±0.071 


0.119 


Z(TIR LVL ) 


£ (PAH 8 fim) 


£(160 


fan) 


0.957 ±0.632 


0.040 ±0.028 


0.946 ±0.042 


0.118 


£(™ £08 ) 


£ (PAH 8 fim) 


£(160 


fim) 


0.396 ±1.709 


0.300 ±0.092 


0.713 ±0.122 


0.182 


X(TIR S ub) 


£(24 fim) 


£(160 


fim) 


-0.641 ±0.094 


0.283 ±0.006 


0.753 ±0.008 


0.028 


^(TIRsings) 


£ (24 fim) 


£(160 


fim) 


0.823 ±0.134 


0.435 ±0.019 


0.562 ±0.019 


0.046 


^(TIRlvl) 


£ (24 fim) 


£(160 


fim) 


1.177 ±0.228 


0.452 ±0.017 


0.535 ±0.019 


0.055 


£(77/W 


£ (24 fim) 


£(160 


fim) 


1.955 ±0.602 


0.482 ±0.025 


0.482 ±0.035 


0.067 


Z(TIRsub) 


£ (70 fim) 


£(160 


fim) 


-1.357 ±0.070 


0.325 ±0.006 


0.727 ± 0.007 


0.024 


J] (TIRsings) 


£ (70 fim) 


£(160 


fim) 


0.051 ±0.153 


0.634 ±0.030 


0.377 ± 0.030 


0.051 


^(TIRlvl) 


£ (70 fim) 


£(160 


fim) 


-0.798 ±0.145 


0.571 ±0.013 


0.465 ±0.013 


0.034 


£(77P £08 ) 


£ (70 fim) 


£(160 


fim) 


-1.570±0.635 


0.843 ±0.050 


0.215 ±0.055 


0.075 


Z(TIRsub) 


£(160 fim) 






-3.288 ±0.127 


1.107 ±0.004 




0.049 


£ (TIRsings) 


£(160 fim) 






0.791 ±0.438 


0.989±0.013 




0.152 


T,(TIR LVL ) 


£(160 fim) 






0.457 ±0.496 


0.999 ±0.015 




0.120 


T.(TIR Em ) 


£(160 fim) 






-2.001 ±1.697 


1.081 ±0.050 




0.203 



NOTE. — Coefficient for the fit logy = a + b X log*; + c X log*2- a is the standard deviation of the data 
points around the best fit line. The subscripts refer to the sample. 




8.4 



8.6 
12+log O/H 



9.0 



1.2 r 



1.0 - 



0.8 - 



0.6 
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8.4 8.6 

12+log O/H 



9.0 



FIG. 3. — Top: y-intercept a (left) and slope b (right) of the linear best fit: logL(TIR) = logL(24 fim) + a + b X log (L(8 fim)/L(24 fim)) versus the oxygen 
abundance. The data points are calculated dividing the sample in 5 bins each containing the same number of data points. Here only subregions are taken into 
account. 



creases. This could be because of PAH destruction or the in- 
hibition of PAH emission at 8 fim in regions with intense or 
hard radiation fields. The ratio L(PAH &{im)/L(24 fjm) may 
drop off more gradually in high metallicity galaxies because 



of increased dust extinction. The photons that most strongly 
affect PAHs in those systems can only travel a short distance 
in the ISM, so the drop off in L(PAH 8fim)/L(24 fjm) with 
increasing radiation field intensity is more gradual. In low 
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TABLE 3 

TIR ESTIMATIONS FROM LUMINOSITIES PER UNIT AREA AND THE OXYGEN ABUNDANCE 



T,(TIR SU b) 

^(TIRsings) 
E(77fi £08 ) 

X(TIR S ub) 
£(77/Ws) 
£(77fi £08 ) 

X(TIRsub) 
£ (TIRsings) 
E(77fi £08 ) 

^(TIRsub) 
£(77/Ws) 

X(TIR SUB ) 
J] (TIRsings) 

£(77fi £08 ) 

£(77/W 
£ (TIRsings) 

£(77fi £08 ) 



£ (PAH 8/im) 
£ (PA// 8ftm) 
£ (PA// 8/im) 
£ (PA// 8ftm) 
£ (PA// 8/im) 
£ (PAH 8/im) 

£ (24 jim) 

E (24 fwn) 

S (24 fim) 
E(PA// 
E(PA// 
E(PA// 

E(24 

E(24 

E(24 

E(70 

E(70 

E(70 



E (24 fim) 
E (24 fim) 
E (24 jLtm) 



8 fim) 

8 (im) 

8 (im) 

fim) 

firri) 

fim) 

firri) 

firri) 

fim) 



E (70 fim) 
E (70 fim) 
E (70 /Ltw) 
E (70 /Ltm) 
E (70 fim) 
E (70 /ion) 



5.667±0.138 
3.996 ±0.909 
7.214 ±0.790 
6.122±0.189 
6.283 ±0.828 
8.893 ±1.271 
5.711 ±0.186 
-0.147 ±1.050 
7.522 ±0.982 
3.606 ±0.075 
0.045 ±0.541 
0.690 ±0.461 
3.497 ±0.138 
-1.941 ±0.505 
0.834 ±0.684 
2.450 ±0.161 
-2.037 ±0.491 
-1.042 ±0.682 



0.431 ±0.015 
0.637 ±0.081 
0.336 ±0.066 
0.854 ±0.007 
0.972±0.013 
0.882 ±0.041 
0.805 ±0.007 
0.970 ±0.017 
0.777 ±0.027 
0.393 ±0.006 
0.266 ±0.046 
0.241 ±0.028 
0.326±0.014 
0.054 ±0.065 
0.215 ±0.045 
0.827 ±0.006 
1.004 ±0.008 
1.021 ±0.019 



0.426 ±0.014 
0.339 ±0.081 
0.508 ±0.057 



0.488 ±0.006 
0.732 ±0.048 
0.773 ±0.031 
0.513±0.014 
0.949 ±0.067 
0.755 ±0.058 



0.009 ±0.022 
-0.256 ±0.1 12 
-0.134 ±0.071 
-0.051 ±0.030 
-0.5 14 ±0.106 
-0.471 ±0.099 
0.222 ±0.028 
0.260 ±0.127 
0.092 ± 0.068 
0.122±0.011 
0.069 ± 0.060 
-0.082 ±0.030 
0.308 ±0.018 
0.268 ±0.059 
0.078 ±0.031 
0.453 ±0.022 
0.271 ±0.059 
0.079 ±0.038 



0.062 
0.127 
0.118 
0.086 
0.146 
0.197 
0.084 
0.185 
0.148 



0.031 
0.063 
0.050 
0.056 
0.086 
0.067 
0.069 
0.086 
0.083 



^(TIRsub) 


E (PAH 8 fim) 


£(160 


fim) 


-1.566 ±0.170 


0.164±0.013 


0.939 ±0.017 


-0.178 ±0.015 


0.043 


J] (TIRsings) 


E (PAH 8 fim) 


£(160 


fim) 


4.522 ±0.537 


0.476 ±0.052 


0.5 19 ±0.054 


-0.430 ±0.065 


0.088 


£(77fi £08 ) 


E (PAH 8 fim) 


E(160 


fim) 


2.514 ± 1.504 


0.449 ± 0.084 


0.590±0.105 


-0.357 ±0.079 


0. 1 5 1 


X(TIR S ub) 


E(24|«i) 


£(160 


fim) 


-0.626 ±0.088 


0.274 ±0.006 


0.790 ±0.008 


-0.1 15 ±0.009 


0.026 


X(TIRs,ngs) 


E (24 fim) 


£(160 


fim) 


1.463 ±0.250 


0.419 ±0.018 


0.582 ±0.019 


-0.094 ± 0.032 


0.043 


£(77fi £08 ) 


E (24 fim) 


£(160 


fim) 


2.020 ±0.590 


0.472 ±0.025 


0.504 ±0.037 


-0.055 ±0.032 


0.065 


£(77/W 


E (70 fim) 


£(160 


fim) 


-1.350 ±0.070 


0.328 ±0.006 


0.721 ±0.008 


0.011 ±0.009 


0.024 


E (TIRsings) 


E (70 fim) 


£(160 


fim) 


-0.385 ±0.329 


0.658 ±0.034 


0.351 ±0.034 


0.060 ± 0.040 


0.050 


£(77R £08 ) 


E (70 fim) 


£(160 


fim) 


-1.632 ±0.645 


0.856 ±0.054 


0.198 ±0.061 


0.026 ±0.038 


0.075 


T,(TIR S ub) 


E(160 fim) 






-3. 199 ±0.1 17 


1.151 ±0.005 




-0.180±0.013 


0.045 


X(TIRs,ngs) 


S(160 fim) 






2.848 ±0.792 


1.005 ±0.013 




-0.301 ±0.099 


0.141 


£(77fi £08 ) 


£(160 /^n) 






-1.507 ±1.659 


1.113±0.051 




-0.189 ±0.092 


0.194 



NOTE. — Coefficient for the fit log y = a + b X log.ti +c X log .12 +d X (12 + log O/ZZ) . a is the standard deviation of the data 
points around the best fit line. The subscripts refer to the sample. 



metallicity systems, however, the photons that affect PAH 
emission can propagate much further through the ISM, so 
L(PAH 8/i»i)/L(24 /j,m) decreases much more quickly as the 
intensity of the radiation field (and L(24 /j,m)/L(TIR)) in- 
creases. 

4.1.2. Star formation intensity effect 

As mentioned in section 14.1.11 the metallicity cannot ex- 
plain that some low-metallicity galaxies follow the trend set 
by rather high-metallicity subregions. This hints that another 
parameter is playing a role in driving the correlation. 

In Figure ||] we plot the log {L(PAH 8/xm)/L(24 /urn)) ratio 
color-coded by Y,(TIR), the total infrared luminosity per unit 
area, which is a proxy to the star formation intensity at high 
metallicities. In Figure [5] the fit parameters for the subregions 
are plotted versus I] (TIR). 

We see that the two branches in Figure [4] are separated 
by T,(TIR). Indeed, most galaxies and galaxy subregions 
that have a low Y^(TIR) tend to be on the lower branch, 
the L(24 fim)/L(TIR) ratio is mostly independent of the 
L(PAH 8/im)/L(24 \xrri) ratio. Conversely, whole galaxies and 
galaxy subregions that have a higher Y.(TIR) are on the upper 
branch. The L(24 [im)/L(TIR) ratio becomes dependent on 
the L{PAH 8 /im) /L(24 /urn) ratio on the upper branch com- 
pared to the lower one. That is, for higher T^(TIR), L{24) 
increases as the infrared spectral energy distribution is hot- 
ter. This dichotomy could be due to the "loss" of infrared in 
low metallicity objects. However, most importantly we note 
that for galaxies that have a L(PAH 8/jm)/L(24 /im) ratio typ- 
ically lower than -0.5 dex, the L(24 fim)/L(TIR) ratio de- 
pends almost exclusively on Y^(TIR). Galaxies that have a 
L(PAH 8/ira)/L(24 jim) ratio typically lower than -0.5 dex 




-0.5 0.0 

log(L(PAH 8um)/L(24)) 

FIG. 4. — The axis and the shapes of the symbols are the same as in Figureff] 
The symbols are color-coded as a function of log£(77P) (W kpc~ 2 ). 

happen to have a low metallicity and as such are expected to 
be predominantly on the lower branch. However, the weak- 
ness or even the lack of PAH emission in these galaxies 
implies that the L(24 fim)/ L{TIR) ratio is very sensitive to 
Y,{TIR) at low metallicities. 

Figure [5] made using the same method as for Figure [3] 
confirms quantitatively the results described above. At low 
Y,{TIR) the slope is shallow but increases to reach a lo- 
cal maximum around log £(77/?) ~ 34.3 before decreasing 
slightly. We note that unsurprisingly the behavior is globally 
similar to the one determined for the oxygen abundance. In- 
deed, lower metallicity galaxies are more transparent as they 
contain less dust and therefore the oxygen abundance and 



7 



E (TIR) are partly correlated. 

The o bserved trend also a ccounts for the relation deter- 
mined bv lCalzetti et al.l (12005b (yellow line in Figure [U since 
those authors have fitted star-forming regions in M5 1 remov- 
ing the contribution to TIR from the dust heating evolved pop- 
ulations, maximizing the effect of the star formation intensity. 

4.1.3. Quality of the correlations 

An important point to note is that the scatter around the 
best fit tends to be larger for whole galaxies compared to sub- 
regions. This is likely due to the diversity in the samples. It 
clearly shows that those sa mples present an important variety 
of infrared properties (e.g. Dale et al. 2009). 

We also estimate the TIR luminosity per unit area from the 
24 fim ban d, that has been sh own to be an accurate SFR es- 
timator by Rie ke et ail d2009l) . This is particularly useful in 
case no other mid- or far-infrared wavelength is available. The 
scatter around the best fit is 0.087 dex, yielding an estimate of 
the TIR luminosity with an accuracy slightly over 20%. Un- 
surprisingly, both luminosities per unit area are an increasing 
function of the oxygen abundance (Pearson correlation coef- 
ficient r = 0.45 for the 24 fim emission and r = 0.53 for the 
total infrared emission). The slop e of the correlation - 0.842, 
shallower than what was found bv lRieke et al.l ( 120091) - shows 
that the 24 fim by itself is not a linear tracer of the TIR lumi- 
nosity. Indeed, at higher luminosities, an increasing fraction 
of the infrared emission can be accounted for by the hot dust 
traced by the 24 /im band. 

4.1.4. Summary 

It appears that both the metallicity and Y,(TIR) play a role 
in the determination of the TIR emission from the 8 /im and 
24 /im bands. While those two parameters are intertwined 
(a low metallicity galaxy tends to be more transparent and 
hence have a lower TIR luminosity), they are the dominant 
parameters in different regimes. 

For L(PAH 8/xm)/L(24 fan) < -0.5 dex, 
L(PAH %fan)/L{2A fim) depends only on T,(TIR). 

For L(PAH 8/im)/L(24 fim) > -0.5 dex two branches can 
be seen. For T,(TIR) < 34.3 dex, L(24 fim)/L{TIR) is 
very weakly dependent on L(PAH 8fim)/L(24 fim), higher 
metallicity galaxies and subregions tend to have a higher 
L(PAH 8/im)/L(24 fim) ratio this being due to the increas- 
ingly stronger PAH emission with the metallicity up to a ratio 
of ~ 0.7 dex. For higher metallicity galaxies and subregions, 
an upper branch is described, the L(24 fim)/L{TIR) ratio in- 
creasing with T,(TIR) while the L(PAH 8fim)/L(24 fim) ratio 
decreases, probably due to a combination of both the heating 
of the very small grains which strongly increases the 24 fim 
luminosity and possible a destruction of PAH in intense ra- 
diation fields. In other words, the main parameters are the 
metallicity which drives the strength of the PAH depending 
on the presence of dust carriers, and the star formation inten- 
sity which drives the temperature of the dust controlling the 
luminosity at 24 fim. 

4.2. TIR estimates that include the 70 fim band 

Here we demonstrate how to estimate the TIR luminosity 
combining the 70 fim and 24 fim bands. The fits are presented 
in Figure[6]and the numerical relations in Tables[T]|2]and[3] 

We see that L(70 fim)/L(TIR) depends significantly on 
the metallicity: low metallicity regions tend to have a high 
L(70 fim)/L(TIR) ratio while higher metallicity regions tend 



to have a lower L(70 fim)/L(TIR) ratio. Indeed, low metallic- 
ity galaxies tend to have a bluer L(70 fim)/L(l6Q fim) color. 
This may demonstrate that the color temperature of the "cool" 
(<30 K) dust in creases as met a llicity decreases. This ef fect 
has by shown bv lHeloul (fl986h : lEngelbracht etall (120081) for 
instance. This could occur because the lower metallicity re- 
sults in lower dust extinction, so a greater fraction of the dust 
that is present in low-metallicity systems is heated by strong 
radiation fields. In contrast, dust extinction in high metallic- 
ity systems would tend to make the radiation field appear less 
intense through most of the ISM. Consequently, most of the 
dust in high metallicity systems would be heated by a weaker 
radiation field and would appear cooler than the dust in low 
metallicity systems. The dependence of the fit parameters on 
£(777?) is weak. A similar result has been found by Calzetti 
et al. (2010, submitted). 

In the relation between L(24 fim)/L(70 fim) and 
L(70 fim)/L(TIR), the integrated galaxy colors exhibit a shal- 
lower slope. Interestingly, about 31% of LVL and E08 galax- 
ies have L(70 fim)/L(TIR) values greater than -0.3 dex but 
only 12% of the SINGS galaxies and a mere 0.1% of the 
galaxy subregions have L(70 fim)/L(TIR) values that are this 
high. 

The TIR luminosity per unit area is well correlated with 
the 70 fim one, and the scatter around the best fit is slightly 
smaller than when using the 24 fim band only (0.082 vs 
0.087). This is due to the fact that large dust grains ap- 
pear warmer in low metallicity galaxies. The slope - 0.902 
- slightly under 1 makes the 70 fim is sub-linear tracer of the 
TIR luminosity which traces the star formation rate for galaxy 
subregions but nearly linear for galaxy samples. See Calzetti 
et al. (2010, submitted) and Li et al. (2010, in preparation) 
for the use of the 70 fim band as a star formation tracer. Tak- 
ing into account the oxygen abundance in Table [3] improves 
the precision of the relation. Indeed, as stated earlier, the 
dust temperature is a function of the metallicity. Providing 
the oxygen abundance permits to indirectly take into account 
the dust temperature and therefore correct the estimation of 
the total infrared emission from the 70 fim band only. 

4.3. TIR estimates that include the 160 fim band 

The 160 fim band probes the cold dust, which makes 
it a good proxy to the TIR luminosity in galaxies 
where it is the dominant source of infrared emission 
(see for example iBendo et al.l 120081) . In Figure UJ we 
show the correlation between log (L(24 fim) /L(l 60 fim)) and 

log (L(160 fim) /L(T 1 R)). In Table [TJ we see that the scatter 
around the best fit is smaller than for the 24 fim and the 70 fim. 

Interestingly, neither the metallicity nor £(777?) have 
a significant effect on the parameters of the fit between 
L(24 fim)/L(160fim) and L(l 60 fim) Jh(T I R) as we can see in 
Figure UJ This proves to be an advantage when the metallicity 
of the galaxy is not known or very uncertain. The 160 fim 
luminosity per unit area is tightly correlated with the TIR lu- 
minosity per unit area. Unsurprisingly the slope is slightly 
superlinear (1.107) and the scatter around the best fit is small 
(0.049 dex) making the 160 fim the best MIPS Spitzer band to 
trace the TIR emission along with the 70 fim band. The scatter 
around the best fit is significantly larger for the SINGS, LVL 
and E08 samples however the slope is almost completely lin- 
ear, especially in the case of the SINGS and LVL samples. 
Once again, this is likely due to the diversity of the popula- 
tions constituting the samples. Taking into account the oxy- 
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33.0 33.5 34.0 34.5 35.0 35.5 33.0 33.5 34.0 34.5 35.0 35.5 

log(E(TIR)) [W kpc 2 ] log(L(TIR)) [W kpc 2 ] 

FIG. 5. — y-intercept a (left) and slope b (right) of the linear best fit: \ogL(TIR) = logL(24 fim)+a + blog (L(8 /im)/L(24 /im)) versus T,(TIR). 




-1.0 -0.8 -0.6 -0.4 -0.2 
log(L(24)/L(70)) 



0.0 




-1.0 -0.8 -0.6 -0.4 -0.2 0.0 
log(L(24)/L(70)) 




FIG. 7. — log (L(160 imi)/IXTIR)) versus log (L(24 fim)/L(l60 fim)) . The symbols and the color-coding of the lines are the same as in Figure[2] The color 
codes the oxygen abundance (left) and £ (TIR) (right). 



gen abundance slightly reduces the scatter when estimating 
the TIR emission from the 160 fim one. As for the 70 fim 
band the reason is that the metallicity gives an indication of 
the possible dust temperature. 



As the 70 urn and 160 fim are major contributors to the 
TIR dDale et al.l 120091 Fig. 11c, f), their temperature should 
be a good indicator of the TIR emission. In Figure [8] 
we show the correlation between log (L(70 ^m)/L( 160/um)) 
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and log (L(\6Qiim)/L{TIR)Y We observe that there is a 
tight non-linear relation between the temperature of the 
dust as traced by the L(70 fim)/L(l 60 fim) ratio. At low 
L(70 pm) /L(l 60 fxm), which indicates that the dust is very 
cold, the bulk of the TIR is accounted for by the 160 /mi. A 
larger fraction of the total infrared emission is accounted for 
by the 70 /mi band as the dust temperature increases. It ap- 
pears that unlike the estimation from the 8 /im and 24 /mi only 
for instance, the parameters of the correlation do not seem to 
depend on the the metallicity or Y,(TIR) in any significant 
way. 

5. DISCUSSION 

The relations derived in this paper are dependent on the 
measurement of the total infrared luminosities and oxygen 
abundances. 

W e have estimated th e TIR luminosity using equation [3] 
from lDraine & Lil (120071). In the co ntext of this study, equa- 
tion 4 from bale & Heloul (120021) yields very similar re- 
sults. For all galaxy subregions, (\og{T IRomine / T IRoaie)) = 
0.01 ± 0.02 as can be seen in the left panel of Fig- 
ure [9] The mean offset is similar for integrated SINGS and 
LVL galaxies data but slightly higher for the E08 sample: 
(\og(TIR Dmim /TIR Dale )) = 0.03 ±0.03. The difference in the 
estimate of the TIR luminosity is small compared to the inter- 
nal scatter within the subregions and the integrated galaxies 
which spans about 1 dex in L(24 fim)/L(TIR). Therefore it 
should not affect the results significantly. 

The method used to calculate the oxygen abundance is an- 
other source of uncertainty. The oxygen abundance estimates 
between the PT05 and KK04 are offset by 0.60 ±0.10 dex. 
The shapes of the slope and y-intercept curves versus the oxy- 
gen abundance are not identical, similarly to what was done 
in Figure [3] There are two reasons to this. First of all, the 
two estimators are not in linear relation. Then, the abundance 
gradients derived by Moustakas et al. (2010, in preparation) 
are different for the two estimators as can be seen in the right 
panel of Figure[9] The consequence is that each subregion has 



a different oxygen abundance offset between the two estima- 
tors. The difference in the fit parameters between the two oxy- 
gen abundance estimators for a given observation provides the 
typical uncertainties due to the abundance estimators. How- 
ever, the precise evaluation of the uncertainties would require 
resolved spectroscopic observations at a ~ 6" resolution. If 
we take the difference of the extrema of the fit parameters 
using the KK04 and PT05 estimators, the typical uncertainty 
on the slope can be evaluated to about ±0.2 and is the main 
uncertainty in this study. Even though the uncertainty is im- 
portant, it is significantly smaller than the range spanned by 
the slope, - 1.2 for the PT05 method and 0.8 for the KK04 
one. 

Even if the relations derived here are distance independent, 
the physical scale encompassed by a pixel ranges from 0.7 kpc 
to 3.6 kpc in galaxy subregions. Even for the closest galaxy, 
each pixel encompasses several star forming regions and the 
ISM of the galaxy. However, for the more distant galaxies 
the mixing of the two components increases with the physical 
scale encompassed. To check whether this effect could induce 
a distance dependent bias in our analysis, we plot in FigurefTOl 
the ratio of the luminosity in each infrared band to the total in- 
frared luminosity for each pixel, as a function of the distance 
of the galaxy. We see that there is no obvious trend with the 
distance. 

6. CONCLUSIONS 

Using data for spatially resolved subregions within 13 face- 
on spiral galaxies as well as integrated luminosities for the 
SINGS, LVL, and E08 galaxies, we have derived new rela- 
tions to estimate the total infrared luminosity based on us- 
ing only one or two Spitzer bands, particularly the 8.0 /im 
and 24 /im bands. Relations incorporating 8.0 /im data 
vary significantly with oxygen abundance and especially with 
S (TIR). However, TIR estimates that do not include 8 /j,m 
data are less dependent on oxygen abundances. In particu- 
lar, the relations between the TIR emission and the 70 //m 
or 160 /im bands are relatively independent of oxygen abun- 
dances compared to the 8 /mi and 24 /mi ones. 



APPENDIX 

RECIPES TO CALCULATE THE TIR LUMINOSITY 

We provide here a few general guidelines to use the relations derived in this paper: 

• If available, the 70 /mi and 160 /im bands should be used to estimate the total infrared emission as they provide the most 
accurate result, 

• If only one band is available, the relations provided in Tables [2] and [3] should be used provided the distance is available and 
the target galaxy is resolved, 

• If only one band is available, and an estimation of the oxygen abundance is available, the relations provided in Table [3] 
should be used, 

• If only the 8 /im and 24 /im bands are available, both should be used to determine the TIR emission, in particular when the 
luminosity per unit area is not available. 

We list in Table|4]the most appropriate formulas to be used to derive the total TIR luminosity of a galaxy, based on our results 
parametrized as a function of available information: IR data (single or multiple bands), metallicity, and expected Y.(TIR). 

We thank the anonymous referee for very useful comments that have helped improving and clarifying the paper. This work has 
been supported by NASA ADP grant NNX07AN90G. MB thanks Alexey Vikhlinin for his technical help during the redaction of 
the manuscript. 
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NOTE. — The low, mid and high ranges are only indicative and are defined relatively 
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